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I. INTRODUCTION 

A soliton is a localised wave that propagates without 
change through a nonlinear medium. Such a localized 
wave forms when the dispersion or diffraction associated 
with the finite size of the wave is balanced by the non- 
linear change of the properties of the medium induced 
by the wave itself. Solitons are universal in nature and 
have been identified in such physical systems, as fluids, 
plasmas, solids, matter waves, and classical field theory. 
Spatial optical solitons - self-trapped light beams - have 
been proposed as building blocks in future ultra-fast all- 
optical devices. Spatial solitons can be used to create 
reconfigurable optical circuits that guide other light sig- 
nals. Circuits with complex functionality and all-optical 
switching or processing can then be achieved through the 
evolution and interaction of one or more solitons [lj . This 
concept has now been verified in several optical materials 
0, 0] and a number of new soliton effects have emerged 
through these studies, such as fusion, fission, and forma- 
tion of bound states. Due to the development of materi- 
als with stronger nonlinearities the optical power needed 
to create such virtual circuits has been reduced to the 
milliwatt and even microwatt level, bringing the concept 
nearer to practical implementation. 

The soliton concept is an integral part of studies of 
the coherent excitations of Bose-Einstein condensates 
(BECs) 4j. Such BECs, inherently have a spatially non- 
local nonlinear response due to the finite range of the 
inter-particle interaction potential. Spatial nonlocality, 
which is already an established concept in plasma physics 
|H El- means that the response of the medium at a par- 
ticular point is not determined solely by the wave inten- 
sity at that point (as in local media), but also depends 
on the wave intensity in its vicinity. The nonlocal nature 
often results from a transport process, such as atom dif- 



fusion || , heat transfer [j| 0] or drift of electric charges 
[lT|. It can also be induced by a long-range molecular 
interaction as in nematic liquid crystals which exhibit 
orientational nonlocal nonlinearity [l^ . IT^ . Nonlocality 
is thus a generic feature of a large number of nonlinear 
systems. It has also recently become important in optics 
[01 lisl ITtH . Although nonlocality can have a consider- 
able impact on many nonlinear phenomena, studies of 
nonlocal nonlinear effects are still in their infancy. 

In this paper we will discuss the role of spatial nonlo- 
cality in phenomena associated with propagation of op- 
tical beams in nonlinear media. These include modula- 
tional instability of plane waves, structural stability of 
localized beams and interaction of solitons. 



II. MODEL 

We consider an optical beam propagating along the 
z-axis of a nonlinear medium with the scalar electric 
field E(f, z) = ip(r, z) exp{iKz — iflt) + c.c. Here r spans 
a D-dimensional transverse coordinate space, K is the 
wavenumber, is the optical frequency, and tp(f, z) is the 
slowly varying amplitude. We assume that the refractive 
index change N(I) induced by the beam with intensity 
I(r,z) = \ip(r,z)\ 2 can be described by the phenomeno- 
logical nonlocal model 



N(I) = s / R(?-r)I(£,z)dH, 



(1) 



where the integral J dr is over all transverse dimen- 
sions and s — 1 (s = —1) corresponds to a focusing 
(defocusing) nonlinearity. The response function R(r), 
which is assumed to be real, localized and symmetric 
(i.e. R(r) — R(r), where r = \f\), satisfies the normal- 
ization condition jR(f)df=l. This model of nonlinear- 
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FIG. 1: Various degrees of nonlocality, as given by the width 
of the nonlocal response function R(x) and the intensity pro- 
file of the beam I(x). Shown are (a) local, (b) weakly nonlo- 
cal, (c) general (nonlocal) and (d) highly nonlocal responses. 



ity leads to the following nonlocal nonlinear Schrodinger 
(NLS) equation governing the evolution of the beam 

id z ip + iviV + N(I)ip = 0. (2) 

The width of the response function R(r) determines 
the degree of nonlocality. For a singular response, R(r) = 
S(r) [see Fig^a)], the refractive index change becomes 
a local function of the light intensity, N(I) — sl(f,z), 
i.e., the refractive index change at a given point is solely 
determined by the light intensity at that very point, and 
Eq.(|2J simplifies to the standard NLS equation 

id z ^ + ^VliP + s ^\ 2 = 0. (3) 

With increasing width of R(r) the light intensity in the 
vicinity of the point r also contributes to the index change 
at that point [see Fig^. While Eq. Q is a phenomeno- 
logical model, it nevertheless adequately describes several 
physical situations where the nonlocal nonlinear response 
is due to various transport effects, such as heat conduc- 
tion or diffusion of molecules or atoms. 

Before proceeding further it is worthwhile to mention 
two important physical situations where the convolution 
term in Eq. (|2J) can be represented in a simplified form, 
which allows for an extensive analytical treatment of the 
resulting equation. First is the so-called weak nonlocality 
limit when the width of the response function is much less 
than the spatial extent of the beam [see Fig^b)]. In this 
case one can formally expand 1(f) in a Taylor series and 
retain only the first significant terms, which gives the 
following simple form of the nonlinearity 

N(I) = s(I + 7 ViZ), 7 = ^ / r 2 R(r)df, (4) 

where 7 is a measure of the strength of the nonlocal- 
ity. Here the nonlocal contribution to the Kerr-type local 



nonlinearity is reflected by the presence of the Laplacian 
of the wave intensity. A nonlinearity in this particular 
form appears in fact naturally in the theory of nonlinear 
effects in plasma 7] . It has been shown recently that the 
one-dimensional version of Eq. with nonlinearity Q 
supports propagation of stable bright and dark solitons 
[l7| . Another limiting case, the so-called highly nonlocal 
limit, refers to the situation when the nonlocal response 
function is much wider than the beam itself [see Fig^d)] . 
In this case 

N(I) = sR(r)P, P = J Idf, (5) 

where P is the total power of the beam. Interestingly, in 
this case the propagation equation becomes local and lin- 
ear. It describes the evolution of an optical beam trapped 
in an effective waveguide structure with the profile given 
by the nonlocal response function. This highly nonlocal 
limit was first explored by Snyder and Mitchell in the 
context of the "accessible solitons" [l4| . 

Although it is quite apparent in several physical sit- 
uations that the nonlinear response in general is non- 
local (as in the case of thermal lensing), the nonlocal 
contribution to the refractive index change is often ne- 
glected 0, 0. This is justified if the spatial scale of 
the beam remains large compared to the characteristic 
response length of the medium (given by the width of 
the response function). However, for very narrow beams 
or beams with fine spatial features (such as dark and 
bright solitons) the nonlocality can be of crucial impor- 
tance and should be taken into account in the ensuing 
theoretical model. 

III. MODULATION AL INSTABILITY 

Modulational instability (MI) constitutes one of the 
most fundamental effects associated with wave propa- 
gation in nonlinear media. It signifies the exponential 
growth of a weak perturbation of the wave as it prop- 
agates. The gain leads to amplification of sidebands, 
which break up the otherwise uniform wave and generate 
fine localized structures (filamentation). Thus, it may 
act as a precursor for the formation of bright solitons. 
Conversely, the existence of stable dark solitons requires 
the absence of MI of the constant intensity background. 
The phenomenon of MI has been identified and studied 
in various physical s yste ms, such as fluids [2(|, plasma 
[2l| , nonlinear optics |22j| , discrete nonlinear syste ms J23| , 
waveguide arrays, and Fermi-resonant interfaces |24|. It 
has been shown that MI is strongly affected by mecha- 
nisms such as saturation of the nonlinearity |25| . coher- 
ence properties of optical beams j2(J , and linear and non- 
linear gratings ^27] . The model J2J) permits plane wave 
solutions of the form 

ip(r, z) = ^o^exp(z/c • r- i/3z), (6) 

where po > is the wave intensity. The parameters sat- 
isfy the same dispersion relation (3 = — spo as for 
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the standard local NLS equation, because the response 
function is normalized. We now consider perturbed plane 
wave solutions in the form 

i>{r,z) = h/Po + ai(<f, z)} exp(ifc ■ f- i/3z), (7) 
where 

ai(^,z) = J 'ai(fc) exp(ifc • £ + \z)dk, £ = r — koz, (8) 

is the complex amplitude of the small perturbation re- 
ferred to a coordinate frame moving with the group ve- 
locity kg. One can show that A satisfies the equation 

M 



-k p 



4p 



sR(k) 



(9) 



where k = \k\ denotes the spatial frequency, R(k) is the 
Fourier spectrum of R(r). The sign of A 2 determines 
the stability of the solution. For A 2 > the pertur- 
bation grows exponentially during propagation with the 
growth rate or gain given by Re{A}, indicating MI. Equa- 
tion © shows that the stability properties of the plane 
wave solutions are completely determined by the spectral 
properties of the nonlocal response function. A detailed 
analysis of all the possible scenarios is given in (2^| . Here 
we summarize the main points. 

The Fourier spectrum of typical response functions, 
such as Gaussians, Lorentzians, and exponentials, is pos- 
itive definite. Therefore, for a defocusing nonlinearity 
(s=— 1), A 2 < and the plane wave solutions are always 
stable. In a focusing medium (s=+l) we always have 
A 2 > 0, in a certain wavenumber band symmetrically 
centered about the origin, where k is sufficiently small. 
This means that the system will always exhibit a long 
wave MI in the focusing case, independent of the precise 
details of the response function. Nonlocality tends to 
suppress the instability, decreasing the growth rate and 
the width of the instability band. However, nonlocality 
can never eliminate the instability completely. This effect 
is illustrated in Fig.|2Ia), where we show the instability 
growth rate in the case of a one-dimensional Gaussian 
response function R(x) = exp(— x 2 /a 2 )/(y/Tra). 

A drastically different behavior is observed for response 
functions whose spectrum is not sign definite. Let us 
consider the generalized Lorentzian 30] 



R(x) 
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(10) 



which for large n approximates the rectangular func- 
tion. The spectrum of this function, which for n=2 
is R(k) = exp(-cr|fc|/V2) [cos(ojfc|/\/2) + sm(a\k\/V2)] 
exhibits periodic changes of its sign. For focusing non- 
linearity this may lead to the appearance of higher 
order instability bands in sufficiently nonlocal media 
[see Fig. Hfb) for er=120]. Even more dramatic is the 
behaviour of the plane wave solution in a defocusing 
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FIG. 2: MI gain profile in self- focusing nonlocal media with 
Gaussian response function for po=l (a) and with generalized 
Lorentzian response function for n=2 and po = L25 (b). 
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FIG. 3: MI in self-defocusing nonlocal media with response 
function 1101 . (a) MI gain for n=2 and po — l. (b) Develop- 
ment of instability for n=6 and po=3. 



medium, where this type of nonlocality may promote a 
high frequency instability of the otherwise stable wave 
[see Fig. EJa)]. In Fig.^Jb) we show the development of 
the instability of the plane wave. The response function 
(fT~0f> has recently been discussed in the context of MI of 
partially coherent waves, where it represents the degree 
of coherence of the plane wave. Interestingly, it also pro- 
moted MI in that case |3(]]. It is worth mentioning that 
the spectrum of llOfl for any value of the exponent n > 2 
contains an oscillatory factor responsible for the forma- 
tion of higher order gain bands in a way analogous to the 
case n = 2. 

The above-mentioned stability properties of the plane 
wave may seem surprising in view of the fact that the 
nonlocality acts to smooth out any sharp modulations 
of the wavefront. This is a generic property of the non- 
locality independent of the particular functional repre- 
sentation. One would therefore naively expect identical 
stability properties for all physically reasonable response 
functions. However, one can also look at the action of the 
nonlocality from a different perspective. In the Fourier 
domain the nonlocality acts as a filter with variable trans- 
mission determined by the form of the spectrum of the 
response function. For many nonlocal models, such as 
that with a Gaussian response function, the transmis- 
sion characteristic of the filter has the form of a well- 
behaving, sign-definite function. However, in cases, such 
as those modelled by (|10|) . this filter not only modulates 
the amplitude of the spectral components of the signal 



beam (perturbation to the plane wave), but also inverts 
the phase of selected components (bands). As this in- 
version is equivalent to a change of the sign of the non- 
linearity (say, from defocusing to focusing), it leads, for 
instance, to amplification of certain harmonics in a de- 
focusing nonlinear media. Recently Peccianti et al. con- 
ducted the first experimental studies of MI in a nonlocal 
medium |3l| . They used nematic liquid crystals, which 
are known to exhibit orientational nonlocal nonlinearity 
with a sign-definite exponential response. The effect of 
suppression of MI was clearly demonstrated. 

IV. BEAM COLLAPSE 
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FIG. 4: Propagation of a vortex beam (charge 1) in a self- 
focusing nonlocal Kerr-like medium. The nonlocality param- 
eter a = (a), a = 1 (b), a = 10, (c). 



Collapse is a phenomenon well-known in the theory of 
wave propagation in nonlinear focusing media. It refers 
to the situation when strong self-focusing of a beam leads 
to a catastrophic increase (blow-up) of its intensity in a 
finite time or after a finite propagation distance 133. 133 . 
|34| . Collapse has been observed in plasma waves |35|l. 
electromagnetic waves or laser beams |36|, BEC's |37| . 
and even capillary-gravity waves on deep water j^l ■ 

Strictly speaking, the existence of collapse is an arte- 
fact of the model equation and signals the limit of its 
applicability. In the vicinity of the collapse regime some 
additional (unaccounted for) physical processes will come 
into play and stop the blow-up. Nevertheless, "collapse- 
like" (or quasi-collapse) dynamics can still occur in real 
physical systems when nonlinearity leads to strong en- 
ergy localisation. In fact, recent experiments with ultra 
cold gases provided clear signatures of the collapse-like 
dynamics of atomic condensates |3T, 39]. 

In the past there have been a few attempts to deter- 
mine the role of nonlocality in development of the col- 
lapse of finite-size beams. Turitsyn was the first to prove 
analytically the arrest of collapse for a specific choice 
of the nonlocal nonlinear response |4jj. Recently Perez- 
Garcia et al. discussed the collapse suppression in the 
case of weak nonlocality, in a direct application for a 
BEC [4Lj . The analysis of the collapse conditions in case 
of general BEC nonlocal model is difficult and has so far 
only been done numerically |42j. 

Here we present an analytical approach to beam col- 
lapse in nonlocal media, which is based on the technique 
introduced in |43j , but which extends and generalizes the 
results to be valid for much more general response func- 
tions. We consider the general case of symmetric, but 
otherwise arbitrarily shaped, non-singular response func- 
tions and prove rigorously that a collapse cannot occur. 

For localized or periodic solutions Eqs. and J5J con- 
serve the power (in optics) or number of atoms (for BEC) 
P and the Hamiltonian H, 

P = J W H = ~\J NIdr, (11) 

where HV'lll = / \ip\ 2 dr. In the local limit when the 
response function is a delta-function, the nonlinear re- 



sponse has the form N(I) — si, as in local optical Kerr 
media described by the conventional NLS equation and in 
BEC's described by the standard Gross-Pitaevskii equa- 
tion. In this local limit multidimensional beams with a 
power higher than a certain critical value will experience 
unbounded self-focusing and collapse after a finite prop- 
agation distance. 

It can be easily shown that in the two extreme lim- 
its of a weakly and highly nonlocal nonlinear response 
the collapse is prevented 0, EH- Here we consider 
the general case of symmetric, but otherwise arbitrar- 
ily shaped, non-singular response functions. Introducing 
the D-dimensional Fourier transform and its inverse, it is 
straightforward to show that for N(I) given by Eq. 0J, 
when s=l the following relations hold [Ij] 



\I(k)\ = | / I{r)e lk - r dr\ < P, 
1 



NIdr 



(27T) 



D 



R{k)\I{k)\ 2 dk. (12) 



For any response functions for which the spectrum R(k) 
is absolute integrable, we then have 



NIdr] <P 2 R , Rq = 



(27T) 



D 



\R(k)\dk, (13) 



and hence we get \H\ > \\Vip\\%/2 - R P 2 /2. This in- 
equality shows that the gradient norm 1 1 VV>| || is bounded 
from above by the conserved quantity 2\H\ + R P 2 for all 
symmetric response functions under the only requirement 
that their spectrum is absolute integrable. It represents a 
rigorous proof that a collapse with the wave- amplitude lo- 
cally going to infinity cannot occur in BEC's or Kerr me- 
dia with a nonlocal nonlinear response for any physically 
reasonable response functions, such as Gaussians, expo- 
nentials, sech-shaped, and even generalized Lorentzian 
response functions of the type (|l(Jfl . 

This result represents a considerable generalization of 
the proof presented in |43l ] , which required the spectrum 
of the response function to be strictly positive definite. 

We emphasize that although the nonlocality prevents 
collapse, it does not remove the collapse-like dynamics of 
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high-power beams. To the contrary, as long as the beam 
width is much larger than the extent of the nonlocality 
(width of the response function), the beam will naturally 
contract as in the local case. However the collapse will 
be arrested when the width becomes comparable with 
that of the response function. Finally, during the sta- 
tionary phase of the beam propagation, its width will be 
comparable with the extent of the nonlocality. 

The stabilising effect of the nonlocality can be further 
illustrated by considering the propagation of a vortex 
beam in a self-focusing medium. Such a beam is charac- 
terized by a bright ring with a helical phase front with 
so-called charge, defined as a closed loop contour integral 
of the wave phase modulo 2tt. A typical example is the 
Gaussian-Laguerre beam 

tp(r) = rexp [-(r/r ) 2 ] exp(z^), (14) 

where r and <fi are the radial and angular coordinates, re- 
spectively. This beam represents a vortex of charge one. 
Beams of such structure have been considered as candi- 
dates for vortex-type solitons in nonlinear self-focusing 
media jj^ - However, it is well known that vortex beams 
cannot form stable stationary structures and disintegrate 
rather quickly when launched in self-focusing local Kerr 
nonlinear media |46|. One can expect to achieve a sta- 
bilization of the vortex beam by utilizing the nonlocal 
character of the nonlinearity. If the extent of the non- 
locality is comparable with the size of the vortex beam 
then the resulting refractive index change will form of 
a broad circular waveguide, which could trap the vor- 
tex beam ensuring its stable propagation. In Fig. (@J 
we illustrate the effect of the nonlocality on propagation 
of a vortex beam in the form of Eq. (|14f> . We numer- 
ically integrated Eq. @ assuming a Gaussian nonlocal 
response R(x) — exp(—x 2 /a 2 )/(y/na). The plots in each 
row show the light intensity distribution at various prop- 
agation distances. The upper row represents local Kerr 
medium (cr=0). It is clearly seen that the vortex beam 
experiences strong focusing and quickly breaks up into 
two fragments. As the nonlocality engages (a = 1, mid- 
dle row) the vortex instability is restrained and break-up 
occurs only after longer propagation. The distance of sta- 
ble propagation significantly increases when nonlocality 
becomes larger (er = 10, bottom row). 

V. NONLOCAL STRUCTURE OF 
PARAMETRIC SOLITONS 

Unlike Kerr solitons, the formation of solitons in 
quadratic nonlinear (or x^) materials does not involve 
a change of the refractive index Thus the underly- 
ing physics of quadratic solitons is often obscured by the 
mathematical model. Only recently Assanto and Stege- 
man used the cascading phase shift and parametric gain 
to give an intuitive interpretation of self-focu sing , defo- 
cusing, and soliton formation in x materials JjJ . Nev- 
ertheless several features of quadratic solitons are still 



without a physical interpretation, such as certain struc- 
tural properties and the formation of bound states. We 
will show below that quadratic solitons can be described 
by nonlocal models. Such models provide simple physical 
explanations of these properties and many more, building 
on a simple waveguide analogy [43, H3, LUJ • 

Consider a fundamental wave (FW) and its second har- 
monic (SH) propagating along the z-direction in a x^ 
crystal under conditions for type I phase-matching. The 
normalized dynamical equations for the slowly varying 
envelopes Ei^ 2 (x,z) are then |£2( 

id z E x + d x d 2 x E x + E\E 2 exp{-i(3z) = 0, 
id z E 2 + d 2 d 2 x E 2 + El exp(i/3z) = 0. (15) 

In the spatial domain d\ « 2d 2 , di j2 > 0, and x rep- 
resents a transverse spatial direction. In the temporal 
domain dj.,2 is arbitrary and x represents time. (3 is 
the normalized phase-mismatch. Physical insight into 
Eqs. (jl5(l may be obtained from the cascading limit, in 
which the phase-mismatch is large, /3 _1 — > 0. Writing 
E 2 =e 2 exp(i/3z) and assuming slow variation of e 2 (x,z) 
gives the NLS equation id^x+d^Ex+p-^E^Ei = 0, 
with e 2 —E 2 /(3. However, this model wrongly predicts 
several features that are known not to exist in Eqs. i|15|) 
and even for stationary solutions it is inaccurate, since 
the term d 2 E 2 is neglected 51]. To obtain a more ac- 
curate model we assume a slow variation of the SH field 
e 2 {x, z) in the propagation direction only (i.e., only d z e 2 
is neglected) , which leads to the nonlocal equation for the 
FW 

id z Ei + did 2 x Ei + p^NiEDEl = 0, 

/OG 
R(x-d)E 2 (t;,z)dZ, (16) 
-oo 

with E 2 =(3~ 1 N exp(«/3z). Equations 1)16|) show that the 
interaction between the FW and SH is equivalent to the 
FW propagating in a medium with a nonlocal nonlin- 
earity. In the Fourier domain the response function is a 
Lorentzian R(k)=l/(1 + sa 2 k 2 ), where a=\d 2 / (3] 1 / 2 rep- 
resents the degree of nonlocality and s=sign(d 2 (3). Both 
Eqs. (|15fl and i|16|) are trivially extended to more trans- 
verse dimensions. 

For s=+l, where the x' 2 '-system 1)1 5|l has a family 
of bright (for di>0) and dark (for <ii<0) soliton solu- 
tions R(k) is positive definite and localized, giving 
R(x) = (2cr) _1 exp(— |x|/cr). It is possible to show, e.g., 
that the nonlocal model 116(1 does not allow collapse in 
any physical dimension |5ll ] , a known property of the x^ 2 ) 
system <|15[) not captured by the cascading limit NLS 
equation. The cascading limit fi~ x — > is now seen to 
correspond to the local limit a — > 0, in which the re- 
sponse function becomes a delta function, R(x) — > S(x). 
With the nonlocal analogy one can further assign simple 
physically intuitive models to the weakly nonlocal limit 
(7<Cl and the strongly nonlocal limit cr3>l. 

For s——l R(k) has poles on the real axis and the 
response function becomes oscillatory with the Cauchy 
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principal value R(x) = (2a)~ 1 s'm(\x\/<r). In this case the 
propagation of solitons has a close analogy with the evo- 
lution of a particle in a nonlinear oscillatory potential. In 
fact, it is possible to show that the oscillatory response 
function explains the fact that dark and bright quadratic 
solitons radiate linear waves 

Equations show the important novel result, that 
in contrast to the conventional nonlocal NLS equation 
treated in detail in this work, the nonlocal response of 
the x system depends on the square of the FW, not 
its intensity. Thus, the phase of the FW enters into the 
picture and one cannot directly transfer the known dy- 
namical properties of plane waves and solitons, such as 
stability. The general model (|16|l and its weakly and 
strongly nonlocal limits thus represents novel equations, 
whose properties potentially allow to understand yet un- 
explained dynamical properties of quadratic solitons. 

In contrast, the stationary properties of nonlocal soli- 
tons, such as the dependence of their profiles on ma- 
terial parameters, directly apply to quadratic solitons. 
One can show that the structure of the nonlocal equa- 
tion governing the stationary fields resulting from H16(l 
is identical to the conventional nonlocal model for sta- 
tionary solitons and thus it has the same weakly and 
strongly nonlocal limits with the same exact bright and 
dark soliton solutions. We recently showed that the non- 
local model elegantly explains the structural properties 
of both bright and dark solitons and their bound states 
and that it provides good approximate quadratic soliton 
solutions in large regimes of the parameter space plj . 



VI. INTERACTION OF DARK NONLOCAL 
SOLITONS 

To this point, we have concentrated on the properties 
of individual beams in nonlocal nonlinear media. Because 
of the specific nature of the nonlocality, which results in 
spatial advancing of the nonlinearity far beyond the ac- 
tual spatial location of the beam, it is natural to expect 
strong influence of the nonlocality on interaction of well 
separated localised waves and solitons. For instance, in 
case of two nearby optical beams, each of them will in- 
duce a refractive index change extending into the region 
of the other one, thereby affecting its trajectory. One can 
show that in a self-focusing medium nonlocality always 
provides an attractive force between interacting bright 
solitons. This effect has been recently demonstrated for 
the interaction of bright solitons formed in a liquid crys- 
tal [5J| . It has been shown that even out-of-phase bright 
spatial solitons, which in a local medium always repel, ex- 
perienced strong attraction that can only be overcome by 
a sufficiently large initial divergence of the soliton trajec- 
tories [H^l • As a consequence of the nonlocality induced 
attraction, bound states of out-of-phase bright solitons 
can be formed [H^ . In this section we will describe a 
novel phenomenon of attraction of dark solitons |5fi| in 
nonlocal nonlinear media with a self-defocusing nonlin- 




-10 10 -10 10 
Transverse coordinate 



FIG. 5: Intensity profile of numerically found stationary dark 
nonlocal solitons (solid lines) . (a) Single dark soliton and (b) a 
bound state of two dark solitons for A = 1 and a = 4. Dashed 
lines indicate the soliton-induced refractive index change. 



earity |57j . 

We will concentrate on the interaction of one- 
dimensional solitons. Without loss of generality, we 
consider the exponential response function R(x) = 
(2a-) -1 exp(—\x\/cr). Following the discussion in the pro- 
ceeding section about the equivalence between stationary 
nonlocal and parametric solitons, and the results of Refs. 
|5lL l53j , the nonlocal NLS equation with such a response 
function predicts the existence of single fundamental dark 
soliton solutions with nonmonotonic tails above a certain 
critical value of a. As a result bound states involving two 
or more solitons can be formed. 

In Fig. we show examples of numerically found dark 
soliton solutions and their bound states. Dashed line il- 
lustrates the soliton-induced waveguide structure which 
guides the soliton. These solitons appear to be very ro- 
bust. 

The ability for dark solitons to form bound states 
and their subsequent stability is a direct consequence of 
the nonlocality-induced long range attraction of solitons. 
This effect can be qualitatively explained using the self- 
guiding concept. In a local defocusing medium the re- 
fractive index change corresponding to two distant dark 
solitons has the form of two waveguides separated by a 
region of lower refractive index (a potential barrier). In 
the presence of nonlocality the effect of the convolution 
term in Eq. is to decrease the index difference be- 
tween these two separate waveguides (lower the barrier) 
thereby allowing light to penetrate the area between soli- 
tons. This, consequently, manifests as soliton attraction. 
The attraction of solitons can be clearly demonstrated by 
simulating numerically the interaction dynamics of two 
nearby solitons. The results of these simulations are sum- 
marized in Fig.[|Ja-c). They clearly demonstrate that as 
the nonlocality parameter a becomes comparable with 
the separation of the solitons they strongly attract and 
trap each other and subsequently propagate together as 
a bound state exhibiting transverse oscillations. 
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FIG. 6: Attraction of dark nonlocal solitons formed either by 
two closely spaced phase jumps (a,b) or a dark notch (gap) 
in the initial cw background intensity (c). In (a) the phase 
jump is 7T and the degree of nonlocality is <j=2, while the 
initial soliton separation is xo=5.5, 4, 2.5 from left to right. 
In (b) the phase jump is 0.957T and £0=2.5, while o=0.1, 1, 
2. In (c) the width of the intensity gap is 7.5, while o=0.1, 3, 
6. 

VII. CONCLUSIONS 

In conclusion, this paper has discussed the properties 
of optical beams propagating in nonlocal nonlinear Kerr- 
like media. We have shown that nonlocality leads to a 
variety of novel phenomena. In particular, it modifies 
the stability of plane waves, depending on whether the 
spectrum of the nonlocal response function is positive 
definite or not. If the spectrum is positive definite then 
nonlocality always tends to suppress MI, but can never 
eliminate it completely. If the spectrum is not positive 



definite then nonlocality may lead to the appearance of 
higher order instability bands for focusing media. For 
defocusing media the effect of nonlocality can be even 
more dramatic as it may actually initiate instability of 
the plane waves, leading to the formation of localised 
structures. 

In local Kerr media multi-dimensional beams are un- 
stable and collapse if they have a power above a certain 
critical value. We have shown that a nonlocal nonlincar- 
ity always will arrest a collapse and allow for the forma- 
tion of stable multi-dimensional solitons. Furthermore, 
we have demonstrated that quadratic solitons formed in 
the process of second harmonic generation are equivalent 
to solitons of a nonlocal medium with an exponential re- 
sponse function. Finally, we have shown that nonlocality 
induces attraction of normally repelling dark solitons and 
allows for the formation of stable bound states. 

In this review we have specifically concentrated on sys- 
tems with spatially nonlocal, Kerr responses. However, 
in many physical systems (BECs 0] as well as pulse 
propagation in non- isotropic media) , the Kerr response is 
composed of a nonlocal (delayed) response supplemented 
by a local contribution within a certain ratio. In the 
modelling of ultrashort pulse propagation nonlocality ap- 
plies to the temporal response. In that case the causality 
principle requires the nonlocality to be represented by 
non-symmetric response functions, and hence the results 
presented in this review are not applicable to this situa- 
tion. For more details on temporal aspects of nonlocality 
see [58|, EH EH an d the references therein. 
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